home *** CD-ROM | disk | FTP | other *** search
-
-
-
- ____IIIIMMMMTTTTQQQQLLLL2222((((3333FFFF)))) ____IIIIMMMMTTTTQQQQLLLL2222((((3333FFFF))))
-
-
-
- NNNNAAAAMMMMEEEE
- IMTQL2, SIMTQL2 - EISPACK routine. This subroutine finds the
- eigenvalues and eigenvectors of a SYMMETRIC TRIDIAGONAL matrix by the
- implicit QL method. The eigenvectors of a FULL SYMMETRIC matrix can also
- be found if TRED2 has been used to reduce this full matrix to
- tridiagonal form.
-
-
- SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
- ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee iiiimmmmttttqqqqllll2222((((nnnnmmmm,,,, nnnn,,,, dddd,,,, eeee,,,, zzzz,,,, iiiieeeerrrrrrrr))))
- iiiinnnntttteeeeggggeeeerrrr nnnnmmmm,,,, nnnn,,,, iiiieeeerrrrrrrr
- ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn dddd((((nnnn)))),,,, eeee((((nnnn)))),,,, zzzz((((nnnnmmmm,,,,nnnn))))
-
- ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ssssiiiimmmmttttqqqqllll2222((((nnnnmmmm,,,, nnnn,,,, dddd,,,, eeee,,,, zzzz,,,, iiiieeeerrrrrrrr))))
- iiiinnnntttteeeeggggeeeerrrr nnnnmmmm,,,, nnnn,,,, iiiieeeerrrrrrrr
- rrrreeeeaaaallll dddd((((nnnn)))),,,, eeee((((nnnn)))),,,, zzzz((((nnnnmmmm,,,,nnnn))))
-
-
-
- DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
- On INPUT
-
- NNNNMMMM must be set to the row dimension of two-dimensional array parameters
- as declared in the calling program dimension statement.
-
- NNNN is the order of the matrix.
-
- DDDD contains the diagonal elements of the input matrix.
-
- EEEE contains the subdiagonal elements of the input matrix in its last N-1
- positions. E(1) is arbitrary.
-
- ZZZZ contains the transformation matrix produced in the reduction by TRED2,
- if performed. If the eigenvectors of the tridiagonal matrix are desired,
- Z must contain the identity matrix. On OUTPUT
-
- DDDD contains the eigenvalues in ASCENDING order. If an error exit is made,
- the eigenvalues are correct but UNORDERED for indices 1,2,...,IERR-1.
-
- EEEE has been destroyed.
-
- ZZZZ contains orthonormal eigenvectors of the symmetric tridiagonal (or
- full) matrix. If an error exit is made, Z contains the eigenvectors
- associated with the stored eigenvalues.
-
- IIIIEEEERRRRRRRR is set to ZERO for normal return, J if the J-th
- eigenvalue has not been
- determined after 30 iterations. Calls PYTHAG(A,B) for sqrt(A**2 +
- B**2). Questions and comments should be directed to B. S. Garbow,
- APPLIED MATHEMATICS DIVISION, ARGONNE NATIONAL LABORATORY
-
-
-
-
-
- PPPPaaaaggggeeee 1111
-
-
-
-